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Within the framework of nonrelativistic QCD, we make a detailed discussion on the doubly 
heavy baryon production through the e + e~ annihilation channel, e + e~ — > 7/Z — > Hqq/ + Q + Q', 
at a high luminosity e + e~ collider. Here Q^ 1 ' stands for the heavy b or c quark. In addition 
to the channel through the usually considered diquark state (QQ')[ 3 Si]$, contributions from 
the channels through other same important diquark states such as (QQ')l 1 Sq]g have also been 
discussed. Uncertainties for the total cross sections are predicted by taking m c = 1.80 ± 0.30 
GeV and mb = 5.10 ± 0.40 GeV. At a super Z- factory running around the Z° mass and with a 
high luminosity up to C oc 10 34 ~ 10 36 cm _2 s _1 , we estimate that about 1.1 x 10 5 ~ 7 H cc events, 
2.6 x 10°~ 7 Hft c events and 1.2 x 10 4 ~ 6 Sbb events can be generated in one operation year. Such a 
Z-factory, thus, will provide a good platform for studying the doubly heavy baryons in comparable 
to the CERN large hadronic collider. 

PACS numbers: 13.66.Bc,12.38.Bx,14.20.-c 



I. INTRODUCTION 

Theoretically, the production of the doubly heavy 
baryons Eqq/ has been analyzed in Refs. |lTjl5l|. where 
the symbol stands for the heavy b or c quark accord- 
ingly. In particular, a computer program GENXICC for 
simulating the hadronic production of the E cc , S& c and 
Hfes has been completed and upgraded in Refs. TH 13. 
which is written in a PYTHIA-compatible format 16j . 
Throughout the paper, Eqqi, corresponding to E cc , E bc 
or respectively, is a short notation for the baryon 
Eqq' q with the light quark q equals to u or d or s re- 
spectively 0. Experimentally, among the doubly heavy 
baryons, only E cc has been observed by SELEX collab- 
oration fl7H19j . Neither Babar collaboration nor Belle 
collaboration have found the evidence for E cc in related 
experiments [2(1 Hlj , to say nothing of the S& c and E b b 
baryons. We hope the CERN large hadronic collider 
(LHC), due to its high collision energy and high lumi- 
nosity, will change the present situation, and especially, 
improve our understanding on those baryons' hadronic 
production properties. 

Comparing to the pp, ep and 77 collisions, a e + e~ 
collider is helpful and has some advantages to perform 
precise measurements for certain processes. To seek the 
Eqqi events at the LHC is feasible, but its hadronic 
background is much noisy in comparison to a e + e" col- 
lider. In the hadronic production, the baryons are pro- 
duced through scattering or annihilating or fusion of 
two initial partons inside the incident hadrons. In ad- 
dition to the dominant gluon-gluon fusion mechanism, 
one also needs to take the extrinsic or intrinsic heavy 
quark mechanisms into consideration, especially at the 
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1 In the present paper, we will ignore the isospin-brcaking effect, 



small pt regions 0, [l^. Hence, the hadronic produc- 
tion becomes much more complicated due to the intro- 
duction of non-perturbative parton distribution functions 
and intrinsic components of the hadron (even though 
they are universal). While at the e + e~ collider, one 
only needs to consider the e + e~ annihilation channel, 
e + e~ — > j/Z° -4 Sqq/ + Q + Q' , which allows one to 
study Eqqi baryon's own properties. 

If the luminosity of a e + e~ collider is C 



for instance, 



denotes E , or 

cca 



-ecu 



or flees accordingly. 



•Lq34-36 cto -2 s -i anc j j.j. s colliding energy is around the 
Z°-pe&k (it is called as a super Z factory [13]), it will 
raise the production rate up to several orders in compar- 
ison to the previous LEP and Belle and Babar experi- 
ments. This increment has already been observed in the 
doubly heavy meson production due to the Z°-boson res- 
onance effect [23M3TI I . It is thus natural to estimate that 
such a super Z-factory also opens new opportunities for 
studying the Sqq< baryon properties, such as their spec- 
troscopy, their inclusive and exclusive decays, and etc.. 
In the present paper, we will study the semi-inclusive 
production of doubly heavy baryon Sqq< at the super 
Z-factory. 

With the non-rclativistic QCD (NRQCD) frame- 
work [IH, the production of Sqq/ baryon can be fac- 
torized into two steps: The first step is to produce two 
free heavy-quark pairs QQ and Q'Q 1 ', which is perturba- 
tive QCD (pQCD) calculable. This is due to the fact that 
the intermediate 7, gluon or Z° should be hard enough 
to generate a heavy-quark pair. The second step is to 
make the two heavy quarks Q and Q' into a bounding 
diquark (QQ') in [ 3 <Si] (or [ 1 So]) s P m state and in 3 
(or 6) color state accordingly; then it will be hadronized 
into Sqq/ baryon, whose probability is described by the 
NRQCD matrix element. More explicitly, the intermedi- 
ate diquarks in S cc and have two spin-and-color con- 
figurations [ 3 Si]3 and [ 1 <S'o]6; while for the intermediate 
diquark (be) in Sf, c , there are four spin-and-color config- 
urations E bc [ 3 Si]s, E bc [ 3 Si] 6l Ebc^Sols, and E^Soje- 

It has been observed that, for the hadronic production 
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FIG. 1. Typical Feynman diagrams for e + e~ — > (QQ')[n] + 
Q + Q', where Q and Q' stand for the heavy c or b quark, [n] 
represents the spin-and-color quantum number of the inter- 
mediate diquark-state. Two new Feynman diagrams for the 
case of Q 7^ Q' can be obtained by exchanging quark-lines of 
Q and Q' . 



(QQ')[n]+Q + Q')(0 H (n)), (1) 



where the matrix element (O h (n)) is proportional to the 
inclusive transition probability of the perturbative state 
(QQ')[n\ pair into the heavy baryon Sqq/, [n] repre- 
sents the spin-and-color quantum numbers for the inter- 
mediate diquark state. The short-distance cross-section 
da(e + e~ — > (QQ')[n] + Q + Q') takes the following form, 



da (e+e- -> (QQ')[n\ + Q + Q>) = 



Z\M\ 2 d$ s 



4>/(Pi ' P2) 2 - rnj 



channels, the contributions from other spin and color 
configurations of the diquark such as (QQ')l 1 So]e con- 
figuration and etc. can also provide sizable contribu- 
tions in addition to the dominant (QQOPSiJa @, H, . 
We will show that this is also the case for the present 
considered e + e~ annihilation channel, e + e~ — > y/Z — > 
Sqq/ + Q + Q', thus, one needs to take all these states 
into consideration for a sound estimation. 

The remaining parts of the paper are organized as 
follows. In Sec. II, we present the detailed formulation 
for dealing with the process of e + + e~ — > ~f/Z° — > 
Sqq/ + Q + Q'. In Sec. Ill, we give the numerical re- 
sults and uncertainty discussion. Sec. IV is reserved for a 
summary. 



where Ai is the hard scattering amplitude, ^ means 
we need to average over the spin states of electron and 
positron and to sum over the color and spin of all final 
particles. The three-particle phase space is 
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The phase-space can be generated and integrated with 
the help of the FormCalc progra m [3411 , or a combination 
of RAMBOS [H and VEGAS Hlwhich can be found 
in generators GENXICC Q] 



and BCVEGPY 37 



II. CALCULATION TECHNOLOGY 

Typical Feynman diagrams for e + {p2)e~{p\) — > 
(QQ )[n](qi) + Q(<7 2 ) + Q'(?3) are presented in Fig.©. 
According to NRQCD factorization formulae [33|, its dif- 
ferential cross section can be written in the following fac- 
torization form: 

da(e+e- -> Z QQ , + Q + Q>) 



A. M (e+e" -> (QQ')N + Q + Q') 

The hard scattering amplitude Ai for the diquark pro- 
duction can be related with the familiar meson produc- 
tion through a proper correspondence. 

The hard scattering amplitude M. for the production 

channel, e + (p 2 )e-( Pl ) -> (QQ')N(?i) + 0(<&) + Q'(«s), 
can be written as : 



I 

M ((QQ')N) = u«i ( , jf Q gi ) ^/(fcp-i^Q) • • • s/(fci,m Q )riu S4 {q 2 )B(S, si,sr,qi,M QQ ,) 
\Mqq> J 

*U S2 {^^l^j r 'l S f( k l> m Q') ■ ■ ■ S f( k 'K-l' m Q')K v s 3 (q3) xCxQxVx C rr > ■ (2) 



where Sj is the spin state of the outgoing (anti) quark, 
5 is the spin state of the diquark, and Mqqi is the di- 
quark mass. The parameters Ti, • • •, T p , r' 1; ■ • •, TJ. are 
sequential interaction vortexes (the 7-matrix elements 

only) along the corresponding spinor lines, Sf{kf\ tuqio) 
is the fcrmion propagator in between the interaction ver- 
texes. As a special case, when there is only one inter- 
action vertex (p = 1 or k = 1), there is no fermion 
propagator. B(S, Si, S2] <fi, Mqq/) is the wavefunction of 



(QQ')[n]-diquark. Here, C is the color factor of the pro- 
cess, Q is the gluon propagator, V is the photon or Z a 
propagator, and C rr i is the leptonic part which can be 
expressed as 

C rr ' = v r (p 2 )Tu r >(p 1 ) , 

in which r^'> stands for the spin state of the electron 
(positron), r = 7^ for 7-propagator and T ~ 7 M (j — 
sin 2 9 W — ^7 5 ) for Z°-propagator. 
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As a comparison, the hard scattering amplitude for e + (p 2 )e (p\) —> (Q'Q)[n](qi) + Q{q2) + Q'fe), can be 
the meson (Q'Q)[n] production through the channel, written as : 



M {{Q'Q)[n]) = u Si (q 2 )r 1 Sf{~k 1 ,m Q ) ■ ■ ■ s f (-k p -i, m Q )r p v sl ( ™ Q _ qi ) B(S, si, s 2 ; qi, M qi q) 

\ M Q'Q J 



XU, 



m Q' 
Mnu 



9i I r' 1 s/(fci,m Q /) • ■■s f (k' K _ 1 )T' Ki v S3 (q 3 ) x C' x Q x V x C rr > 



(3) 



where Mq,q is the meson mass, B(S, Si, s 2 ; q%, Mq,q) is A comparison of Eq.© and Eq. ((3J) tells us that if we 
the wavefunction of (<2'Q)[n]-meson an d C is the color can transform the fermion line of the diquark-production 
factor for the present case. Other parameters have the case to an anti-fermion line of the meson-production case, 



same meaning as above. 



then we can derive the diquark amplitude following the 
same way of meson production. More definitely, we need 
to deal with the following fermion line, 



2s i 17 9i T pSf(k p -i,m Q ) ■ ■ ■ Sf(ki,mQ)T 1 v Si (q 2 ). 

\ M QQ> 



Setting T to be the operation of matrix transportation 
and C = — 17 2 7° to be the charge conjugation matrix, 
we obtain the following equations 

= -u s M, c-rjc = -r h cc- = 1, 



C sJ(ki,m Q )C = Sf(-ki,m,Q), 
c - n r ( m Q ( m Q 



qi = v Sl 



Mqq>J 1 \M qq , 
Then, the fermion line a changes to 



a = a 



v Si (q2Wi s f (ki,m Q ) ■■■s f {k p -i,m Q )T p u Si ^j^H 



viMCC-Ti CC-sj (kx,m Q )CC- ■ ■ ■ CCSj {k p _ u m Q )CC'Y p CC~u Sl 

(-l) (p+1) W S4 (g 2 )riS / (-fcl,TOQ)---S / (-fc p _ 1 ,TO Q )r p W Sl ^Jg-gr^ . 



Mr 



-<h 



Thus, Eq.@ can be transformed as 



M((QQ')[n}) = (-l)^u Si {q 2 )T 1 s f (-k 1 ,m Q ) ■ ■ ■ s f {-k p ^m Q )T p v Sl (-^-qi ) B(S, 81, s 2 ; qi , Mqq.) 

\Mqq< 



(i^^7 ?1 ) ^^l' mQ '^ ' ' ' S /( fc K-l) r «^3(<73 



X C X Q X T> X C r 



(4) 



Comparing Eq.Q with Eq.Q, one can, thus, follow the By taking the Lorentz indexes and the color indexes 
same procedures of meson production to finish the calcu- explicitly, the hard scattering amplitude M. in Eq.((4]) 
lation. 
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can be rewritten as: 



iM{{QQ')[n}) =k^^xD f x L v rr , 



(5) 



5=1 



For the production through Z°-propagator: k 
-C, 



s 2 e w ' 



fc 2 



k^k v 
k 2 



For the production through 7-propagator: the overall pa- with T z the total decay width of Z° boson, and 
rameter K=eQ(/)e 2 g 2 Cy , where eg(/> is the electric charge, _^ 

in unit of e, of the quark and C^- is the color factor C v rr t = v r (p2)^ v ( sin 2 9 W 7 s I u r /(pi). 

with i and j the color indices of the outgoing antiquarks, ^ ' 

the propagator V, v = ^g, v and the leptonic vector ^ vectorg ^ ^ = 1 4) ^ bg read from ^ 

Feynman diagrams in Fig.([T]). More explicitly, Ag can 

C rr i = v r {p2)^ v u r i{pi). be written as 



^1 = u s (q 2 )j p 



Tt ^ In \ 

(911+92) 2 {q! + q 2 ) 2 - m z , z ^ 



^ + $2+ m Q> 



(6) 



n (QQ')[»]( gl ) 



^2 = M s(?2)7p 



-r! 



t 2 - $g 2 - m? Q , 



(7) 



A% = -u s (q 2 )T 



--Q 



(~9i - 93) 2 - m 2 n lp (912 + 93) 



(8) 



and 



g/2 + jf 3 + jf 2 +mq ^ ' 



-^4 = - u s (<72)7p 



(9/2 + ^ + - m 2 ^"(912 + 93) 



(9) 



where for convenience, we introduce a general interaction 
vertex r^ (() , the needed ones are 



r£ c = 7^ 



1^ = 7" 



I 1 2 -2 

a + 8 1 sir 1 

1 4 3 



' I ~ I Sin2 ' 



4 7 



4 r 



the amplitude 0- For the production of (cc)- and (66)- 
diquarks, we only need to calculate and ^2 J but need 
to time the squared amplitude for (cc)- and (66)-diquark 
by an overall factor (2 2 /2!) = 2, where the 1/2! factor is 
due to the symmetry of the diquark wavefunction. While 
for the case of (6c)-diquark production, we need to cal- 
culate A* s (S = l,--,4). 

The projector n?gg,jr, (91), under the non-relativistic 
approximation, takes the following form [33| : 



Here a = 1 and /3 = for 7 — — vertex, and 
a = and /3 = 1 for Z° — QW — QW-vertex, respectively. 
The momenta of the constituent quarks which forms the 
bound state are 



9n 



m Q 
Mqq. 



qi+q and q 12 



M QQ 



■9i - 9, 



(10) 



where Mqq> = tuq + ttiqi is implicitly adopted to en- 
sure the gauge invariance of the hard scattering ampli- 
tude, and q is the relative momentum between the two 
constituent quarks inside the diquark. Due to the non- 
relativistic approximation, q is small and neglected in 



n ?W')M = (£l7 5 + 6/(9l)) (^ + M QQI ), 

where e(9i) is the polarization vector for spin-triplet 
state. Here the projector n 1 is for the case of spin-triplet 
[ 3 5i], which corresponds to £1 =0 and £2 = 1, an d the 
projector H is for the case of spin-singlet I 1 So], which 
corresponds to £1 — 1 and £2 = 0. 



2 The integration over q results in a wavefunction at zero which has 
been absorbed into the overall non-perturbative matrix element. 
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B. The color factor CVj 

Because of the fact that 3(g)3 = 306in SU C (3) color 
group, the (Q<5')-diquark can be either in anti-triplet 3 
or in sextuplet 6 color state. According to Fig. (JTJ) , the 
color factor Cy of the process is 

Cy = N X J2( Ta )im(T a ) jn X G mnk (12) 

where i, j,m,n = 1,2,3 are color indices of the two outgo- 
ing anti-quarks Q and Q' and the two constituent quarks 
Q and Q' of the diquark respectively, k is the color in- 
dices of the diquark (QQ'), and the indices a = 1, • • • , 8 
is the color index for the gluon, and the normalization 
constant Af = \/T/2. The function G mn k equals to the 
anti-symmetric function e m jk and equals to the symmet- 
ric function f m jk accordingly: the anti-symmetric e m jk 
satisfies 

and the symmetric f m jk satisfies 

fmjkfvn'j'k — $mm' $ j ]' $mj' $jm' • 

Then, the square of Cfj equals to | for the color anti- 
triplet diquark production, and | for the color sextuplet 
diquark production, respectively. 

C. Transition from (QQ')[n] to Hqq/ 

According to NRQCD, the Sqq/ can be expanded as a 
series of Fock states which are according to the relative 
velocity (v) of the constituent heavy quarks in the baryon 
rest frame : 

\3QQ>) = Ci\(QQ')q) + c 2 \(QQ')qg) + c 3 \(QQ')qgg) + ■ ■ ■ , 

where Ci(i = 1, 2, ■ ■ •) is a function of the small velocity 
v, where v is the relative velocity of the heavy quarks in 
the rest frame of the diquark. As for S cc and S&;,, the 
intermediate binding diquark in the Fock states can be in 
either (QQ^S^ or {QQ^S^e (Q = c,b) state; while 
for S(, c , the intermediate binding diquark within the Fock 
states can be in (bc)[ 3 Si]s or (&c)P_Sb]a or {bc)[ 3 Si]e or 
(&c)[ 1 5o]6 states respectively. Ref.Q has suggested that 
the baryon can be formed with the component \{QQ')qg) 
as well as the usual \{QQ')q): One of the heavy quarks 
emits a gluon, which does not change the spin of the 
heavy quark, and this gluon splits into a qq, the light 
quarks can also emit gluons, then the component can 
be formed with the light quark q plus one gluon. Be- 
cause a light quark can emit gluons easily, one can take 
the transition probability from these diquark states into 
the corresponding baryon to be at the same importance. 
Then, as a rough order estimation, one can take the tran- 
sition probability from these diquark states into the cor- 
responding baryon to be the same; i.e. 

h 6 ^ hs, (13) 



where h$ stands for the probability of transforming the 
color anti-triplet diquark into the baryon and he stands 
for the probability of transforming the color sextuplet 
diquark into the baryon. These non-perturbative matrix 
elements can be determined from the potential model or 
from the non-perturbative methods such as QCD sum 
rules and lattice QCD. 

It is also noted in Ref. Q that if the baryon is formed 
by the Fock state component \{QQ')q) only, the emit- 
ted gluon for the case of {QQ 1 ) in 1 5'o spin-state, which 
will split into a gg-pair with q being combined by {QQ') 
to form the baryon, must change the spin of the heavy 
quark. Then, according to NRQCD [321, one may con- 
clude that he must be i^-suppressed in comparison with 
h%. This provides the underlying reason why only h<$ 
has been taken into consideration for the doubly heavy 
baryon production, c.f. Refs. [H-Q. Fortunately, these 
matrix elements are overall parameters and their uncer- 
tainties can be conveniently discussed. We will adopt the 
approximation (| 13[) throughout the paper 0. 

Furthermore, the color-singlet nonperturbative matrix 
element (O h (1S)) in the factorization formulae (U} can 
be related to the Schrodingcr wavefunctions at the origin 

i^qq<)(o)i m 

he^h- 3 = (<D H {IS)) ~ |V|(QQ')[is]>(0)| 2 - (14) 

Since the spin-splitting effect is small, we do not dis- 
tinguish the bound-state parameters for the spin-singlet 
and the spin-triplet states; i.e. those parameters, such as 
the constituent quark masses, the bound state mass, the 
wavefunction and etc, are taken to be the same for the 
spin-singlet and the spin-triplet states. 

III. NUMERICAL RESULTS 

When doing the numerical calculation, the input pa- 
rameters are taken as the following values [9| : 

|* cc (0)| 2 = 0.039 GcV 3 , 
|* bc (0)| 2 = 0.065 GcV 3 , 
|* 66 (0)| 2 =0.152 GcV 3 , 
m c = 1.8 GcV, m b = 5.1 GcV, 

M Bcc = 3.6 GeV, M Ebc = 6.9 GeV, M 3bb = 10.2 GeV . 

Other input parameters are taken from the Particle Data 
Group H: r * = 2 - 4952 GeV : m z = 91.1876 GeV, 
mw = 80.399 GeV, and cos8 w = mw/mz- We take 
the renormalization scale to be 2m c for S cc or S& c , 



3 In different to the heavy quark fragmentation, during the frag- 
mentation of diquark into baryon, the diquark may dissociate, 
which will decrease the baryon production cross section to a cer- 
tain degree. In our present calculation we will not take this effect 
into consideration. Thus, our present estimations can be treated 
as an upper limit for the total cross sections. 
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cross section (pb) 


e+e" -> 7 -> S^P-Sils) 


8.90 x 10- 4 


e+e" ->• 7 -»• Seethe) 


4.29 x 10- 4 


e + e - ^ z° -»• H cc ([ 3 Si]a) 


0.727 


e+e" -> 2° -> Scc^SoJa) 


0.353 


e+e" -> 7 -> Hi, c ([ 3 Si]g) 


2.66 x 10- 4 


e+e" -> 7 ->■ H 6c ([ 3 Si] 6 ) 


1.33 x 10" 4 


e+e" ->• 7 -> Sbcd^SoJs) 


1.85 x 10- 4 


e+e" -)• 7 -)• Sbcf^SoJe) 


9.27 x 10" 5 


e +e- Z° ->■ SteCpSJa) 


1.00 


e+ e - ^ Z° ^ Hi, c ([ 3 Si]a) 


0.502 


e+e" -> Z° ->• H^Sols) 


0.730 


e+e" -> Z° -> Enable) 


0.365 


e+e" -> 7 -> S M ([ 3 Si]g) 


1.94 x 10~ 5 


e+e" ->• 7 -> S i)1) ([ 1 5o]e) 


8.95 x 10" 6 


e+ e - -> Z° ->■ SwtpSila) 


8.03 x 10~ 2 


e+e" -»• Z° -S- SwCfSoJa) 


3.86 x 10" 2 



TABLE I. Total cross section (in unit: pb) for the production 
of S,QQ/(n) baryons (n stands for the intermediate diquark 
state) through e+e" annihilation at the Z°-Peak (yS = mz). 



and 2nib for 3^, which leads to a s (2m c ) = 0.212 and 
a s (2m(,) = 0.164 for leading-order a s running. 

At the super Z-factory, total cross sections for the pro- 
duction channel, e+ + e~ — > Eqq> (n) + Q + Q' (n stands 
for the intermediate diquark state) , are presented in Ta- 
ble |H It is found that, 

• When the e+e" collision energy is at the Z°- 
peak {E cm = y/S = mz), the production channel 
through a 7-propagator is much smaller than the 
case of a Z°-propagator; i.e. for the production 
through the same diquark state, its cross-section is 
always less than 10 -3 of the ^°-case. 

• At the Z°-peak, the contribution from the diquark 
state n = pSole configuration is almost half of 
that of n = [ 3 Si]s for the E cc and Ebb baryon 
production. While for the Sf, c baryon production, 
four intermediate (6c)-diquark states will make siz- 
able contributions: the total cross-sections of n = 
[ 3 Si] 6 , P-Sola and pSoJe are about 50%, 73% and 
37% of that of n = [ 3 S±\s. Then, one should take 
all these Fock states' contributions into considera- 
tion in order to derive a sound estimation. 



To show how the total cross sections change with 
the variation of e+e - collision energy E cm = y/S, we 
present their cross sections versus the collision energy 
in Figs. (|2l3l4j) for the E cc , Eb c and Ebb baryon produc- 
tion respectively. In small collision energy region (e.g. 
E cm ^5 50 GeV) , we have 




40 50 60 70 

E GeV 

cm 



FIG. 2. Total cross sections (in unit: pb) of the production 
channel, e+ + e~ — > j/Z° — > H cc (n) + c + c, versus the e+e" 
collision energy E cm = y/S, where n stands for the corre- 
sponding intermediate (cc)-diquark state. 



Ecm, 


95%m z 


97%m z 


m z 


103%m z 


105%m z 


ff H cc ([30i] H ) 


0.0487 


0.1215 


0.7283 


0.1307 


0.0543 


°"H cc ([lSol 6 ) 


0.0236 


0.0589 


0.3533 


0.0634 


0.0263 




0.0658 


0.166 


1.00 


0.180 


0.0744 




0.0329 


0.0830 


0.502 


0.0899 


0.0372 


a H bc ([lS ] s ) 


0.0478 


0.121 


0.730 


0.131 


0.0542 




0.0239 


0.0603 


0.365 


0.0655 


0.0271 




0.0052 


0.0132 


0.0803 


0.0145 


0.0060 


^([iSole) 


0.0025 


0.0063 


0.0386 


0.0070 


0.0029 



3/ (re) > a e+e- 



'CI- 



TABLE II. The cross section (in unit: pb) for the Hqq/ baryon 
through the process e+e" — > "f/Z° — > Hqq/(ti) + Q + Q' 
(n stands for the intermediate diquark state) under several 
typical collision energies. 



The production cross section for the case of 7-propagator 
dominates at small E cm , there is a small peak at 10 ~ 30 
GeV, and then, it decreases slowly with the increment 
of E cm . The production cross section for the case of 
Z°-propagator is negligible at small E cm , but it arises 
logarithmically, and when E cm ~ mz, due to the Z°- 
resonance effect, it is about four-orders bigger than its 
value at E cm = 50 GeV and is about three-orders bigger 
than the total cross-section (summed up cross-section for 
both the 7-propagator and the Z°-propagator) in low 
region of E cm ~ 10 — 30 GeV. This shows that a e+e" 
collidcr running around the Z°-peak will provide us a 
better chance for producing more doubly heavy baryons 
than the previous LEP, Babar and Belle platforms. 

To be useful reference, we calculate the total cross sec- 
tion (summed up cross-section for both the 7-propagator 
and the Z°-propagator) by taking the collision energy 
E cm = (1 ± 3%)mz and E cm = (1 ± 5%)m z , which 
are shown by Table HU It shows that when the collision 
energy E crn is away from its center value of mz, the to- 
tal cross sections drops down quickly: a 5% deviation 
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FIG. 3. Total cross sections (in unit: pb) of the production 
channel, e + + e~ — > j/Z° — > Hb c (n) + c+b, versus the e + e~ 
collision energy E cm = \i r S, where n stands for the corre- 
sponding intermediate (foc)-diquark state. 



will lead to about one-order-lower cross section. Tak- 
ing the production of S cc [ 3 S'i]3 as an example, the total 
cross section for E cm = (1 ± 3%)mz will be lowered to 
~ Qg'Sw) of its peak value; and the total cross section 

for E cm = (1 ± 5%)mz will be lowered to ~ (57%) °f lts 
peak value. 

Considering the super Z-factory which is running 
around the mass of Z°-boson and with a high luminos- 
ity C oc 10 34 cm~ 2 s _1 , the number of Eqq'(Q,Q' = c, b) 
events per year can be estimated : 

• Totally, 1.08 x 10 5 S cc -baryon events can be gen- 
erated, which includes 7.28 x 10 4 events coming 
from the intermediate (cc)[ 3 Si]q diquark state and 
3.53 x 10 4 events coming from the intermediate 
(cc)[ 1 So}e diquark state. 

• Totally, 2.60 x 10 5 5f, c -baryon events can be gener- 
ated, which includes 1.00 x 10 5 events coming from 



4 Approximately, 1 year fa n X 10 7 s, but it is common that a 
collider only operates about I /it of the time a year |39( . i.e. we 
are customary to take 10 34 cm -2 s — 1 fa 10 5 pb — 1 /year. 



FIG. 4. Total cross sections (in unit: pb) of the production 
channel, e + + e~ — > — > Sbb(n) + b + b, versus the e + e~ 
collision energy E cm = \f~S, where n stands for the corre- 
sponding intermediate (&6)-diquark state. 



m c (GeV) 


1.50 


1.65 


1.80 


1.95 


2.10 


°"S co ([3Sl] B ) 


1.270 


0.949 


0.727 


0.569 


0.453 


ff B cc (FiSnl 6 ) 


0.616 


0.461 


0.353 


0.276 


0.220 




1.85 


1.34 


1.00 


0.767 


0.598 


^H hr ([3 Sl ] 6 ) 


0.926 


0.672 


0.502 


0.383 


0.299 


^HbctPSols) 


1.28 


0.954 


0.730 


0.571 


0.455 




0.639 


0.477 


0.365 


0.285 


0.227 



TABLE III. Uncertainties for the total cross section (in unit: 
pb) of the production channel e + e~ — > Z° — > Sqq/ (n)+Q+Q' 
(n stands for the intermediate diquark state) with varying m c . 
Here, mj is fixed to be 5.10 GeV. 



the intermediate (be) [ 3 Si]^ diquark state, 5.02 x 10 4 
events coming from the intermediate (bc)[ 3 Si]e di- 
quark state, 7.30 x 10 4 events coming from the in- 
termediate (be) [ 1 S'o]3 diquark state and 3.65 x 10 4 
events coming from the intermediate (bc)[ 1 So]e di- 
quark state. 

• Totally, 1.19 x 10 4 5f,(,-baryon events can be gen- 
erated, which includes 8.03 x 10 3 events coming 
from the intermediate (bb)[ 3 Si}^ diquark state and 
3.86 x 10 3 events coming from the intermediate 
(bb) I 1 So] 6 diquark state. 

• About 10 4 — 10 5 -order events per year shows clearly 
that sizable events of S cc , 3;, c , and even can be 
produced at the future super Z-factory. If its lu- 
minosity can be increased up to £ oc 10 36 cm _2 s _1 , 
the event numbers will be further increased by two 
orders. Thus, even if we have a slight deviation of 
E cm from mz, e.g. up to 5%-departure, there is 
still sizable events. 

Finally, we discuss the theoretical uncertainties from 
the heavy quark masses by varying m c = 1.80 ± 0.30 
GeV and mb = 5.10 ± 0.40 GeV. The cross sections for 
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m b (GeV) 


4.70 


4.90 


5.10 


5.30 


5.50 


t7 H bo ([3 Sl ] 3 ) 


0.987 


0.995 


1.00 


1.01 


1.02 


^B hr .(\ 3 S l ] e ) 


0.494 


0.498 


0.502 


0.506 


0.509 


^([ISols) 


0.738 


0.734 


0.730 


0.726 


0.723 




0.369 


0.367 


0.365 


0.363 


0.362 




0.105 


0.092 


0.080 


0.071 


0.062 


^([iSole) 


0.051 


0.044 


0.039 


0.034 


0.030 



TABLE IV. Uncertainties for the total cross section (in unit: 
pb) of the production channel e + e~ — > Z° — > EQQi(n)+Q+Q' 
(n stands for the intermediate diquark state) with varying m D . 
Here, m c is fixed to be 1.80 GeV. 



E C m = mz with varying m c and mb are presented in 
Tables IIIII and IIV1 which are more sensitive to m c than 
those of mb- Here we only consider the cross sections for 
the channels through Z°-propagator, since the channels 
through 7-propagator is small and negligible as shown by 
TablclH By adding these two uncertainties in quadrature, 
we obtain 



< 7 e++e-- 


+Z°- 


>a co ([35i] a 


)+cc 


= 0.727+°;|H pb, 


(15) 


C e ++ e -- 


+Z°- 


+E cc ([iSo]e 


)+cc 


= o.353;° pb, 


(16) 


cr e++e-- 


-^Z"- 


>s 6e (ps 1 ] a 


)+bc 


= i.oo;°; f 5 2 pb, 


(17) 


tJ e++e~- 


-+Z - 


^3i,o([ 3 Si] 6 


)+bc 


= 0.502;°;^ pb, 


(18) 


cr e++e-- 


-+z a - 


^H 6c ([iSo]3 


)+bc 


= 0.730;° j 5 75 pb, 


(19) 


a e++e-- 


-+z a - 


^^fac([ 1 5'o]6 


)+bc 


= 0.365;°;^ P b, 


(20) 


a e++e~- 




>H 6i ,([3S 1 ] a 


)+bb 


= 0.080;S;gai Pb, 


(21) 


a e++e-- 


-+z°- 


->E bb ([iS ]6 


)+bb 


= 0.039;°;°°^ pb. 


(22) 



IV. SUMMARY 

Using the NRQCD factorization formula, we have stud- 
ied the production of the doubly heavy baryon Sqq/ (3 cc , 
Sb c , Sbf,) through the e + e~ annihilation. When the e + e~ 
collision energy is around the Z°-peak, sizable cross- 
sections for the doubly heavy baryon can be obtained. 
Contributions from those (<3<3')[n]-diquark states with 
the same importance have been discussed. By adding all 



these intermediate diquark states' contributions together, 
we obtain the total cross section for Sqq/ at E cm = Z° : 





■+7/Z°- 


->H cc +cc 


= 1.08;°^ P b, 


(23) 




-ry/ZO- 


-S-H b<: +bc 


= 2.607.^ P b, 


(24) 


Ce+e-- 


-+~//Z°- 


+5 bb +bb 


= 0.119;°;°^ pb , 


(25) 



where the errors are caused by varying m c = 1.80 ± 0.30 
GeV and m b = 5.10 ± 0.40 GeV. 

At the super Z-factory with a high luminosity up to 
C oc 10 34 — 10 36 cm -2 s -1 , one would expect to accumu- 
late 1.1 x 10 5 ~ 7 3 CC events, 2.6 x 10 5 ~ 7 S 6c events and 
1.2 x 10 4 ~ 6 Ebb events in one operation year. If taking 
the e + e~ collision energy to run slightly off the Z°-peak, 
i.e. \/S = 0.95mz or 1.05m^, the total production cross 
section will be lowered by about one order of magnitude 
from its peak value. However, one may still observe siz- 
able events. At the hadronic colliders, there is much pol- 
lution from the hadronic background and many produced 
baryon events shall be cut off by the trigging condition. 
So, some alternative measurements would be helpful for 
a comprehensive study. In addition to the hadronic col- 
lider LHC, the super Z-factory will provide another good 
platform for studying the Sqq/ baryon properties. 

As a final remark, if not too much baryon events are cut 
off by the trigging condition at the super Z-factory, one 
may have the chance to further distinguish the baryons 
with different light quark content. When an intermediate 
diquark is formed, it will grab a light quark (with gluons 
if necessary) from the 'environment' to form a colorless 
doubly heavy baryon with a relative possibilit y fo r var- 
ious light quarks as u : d : s ~ 1 : 1 : 0.3 [la ]. For 
example, if the diquark (cc) [ 3 5i]g is produced, then it 
will fragment into S+ + with 43% probability, 5+ with 
43% probability and f2+ with 14% probability accord- 
ingly. The condition is the same for the production of 
S° b , 3r/ b and fir/ 6 or the production of 3^, Sg c and f2° c . 
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